One-step replica symmetry breaking solution of the quadrupolar glass model 
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. We consider the quadrupolar glass model with infinite-range random interaction. Introducing a 

' simple one-step replica symmetry breaking ansatz we investigate the para-glass continuous (discon- 

tinuous) transition which occurs below (above) a critical value of the quadrupole dimension m* . By 
q , using a mean-field approximation we study the stability of the one-step replica symmetry breaking 

solution and show that for m > m* there are two transitions. The thermodynamic transition at a 
temperature To is discontinuous but there is no latent heat. At a higher temperature we find the 
dynamical or glass transition temperature To and the corresponding discontinuous jump qG of the 
order parameter. 
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I. INTRODUCTION 



In the last decades quadrupolar glasses have found widespread experimental and theoretical interest . Disordered 
quadrupolar glasses are produced by random dilution of (quadrupolar) molecular crystals with atoms which have no 
quadrupole moments; well-known examples of such systems are K(CN) x Bri_ x or Na(CN) x Ch- x , or N2Ar, CuCN, 
or solid hydrogen (see Ref. 2] for a review). The success of the Sherrington-Kirkpatrick (SK) model 3] in providing a 
good theory to describe systems of interacting magnetic or electric dipole moments, suggests to extend to quadrupolar 
glasses the same kind of analysis. However, there are differences between standard spin glass systems and quadrupolar 
£^ ■ glasses; the latter do not have the global inversion symmetry Si — ► —Si for all spins. For several systems without 
reflection symmetry, and close to the transition temperature, only one step in the Parisi replica symmetry breaking 
scheme is sufficient to describe the transition para-glass above a lower critical dimension Q|. Indeed, the one-step 
replica symmetry breaking (1RSB) scheme has proven to provide stable solutions for the Potts glass model jMI&LjIj 
the spherical p-spin model [8| and the p-spin Ising spin glass model [!| . It is the purpose of the present paper to show 
that the IRSB scheme can be applied also to the quadrupolar glass model, and indeed it provides a stable solution 
in certain regimes. 

In the present paper we shall consider a perturbative evaluation of the free energy by means of a Taylor expansion 
up to fourth order in the order parameter. It is obvious that the perturbative approach is most reliable near the 
transition temperature. We shall show that the transition from the replica symmetric (RS) state to the IRSB occurs 
' either discontinuously or continuously, depending on the value of the quadrupole dimension m. A similar dependence 
of the RS to IRSB transition, though on the value of an external field, is exhibited by the spherical p-spin model || 
' and the Ising p-spin model For any p > 2, the transition is discontinuous (continuous) for fields weaker (larger) 
than a critical value of the external field h c , which depends on p. 

The plan of the paper is as follows: in Sec. II we shall use a rather pedagogical approach mainly to review the 
results obtained in the mean-field analysis of the quadrupolar glass model in the framework of the replica symmetry 
ansatz |To| . Section III is devoted to the study of the IRSB solutions of the saddle-point equations, assuming all the 
transitions to be continuous or at worst weakly discontinuous. In Sec. IV we shall perform the de Almeida-Thouless 
(AT) stability analysis while the dynamical transition is discussed in Sec. V. The concluding remarks are given in 
Sec. VI. 



5— I 1 II. UNIAXIAL QUADRUPOLAR GLASS 

] 

The infinite-ranged quadrupolar glass model has been first introduced by Goldbart and Sherrington (GS) 10]. The 
model assumes the quadrupole-quadrupole interaction to be more dominant than the interactions between dipoles. 
This appears to be the case in several experimental situations, where the quadrupolar species occupy the sites of a 
regular lattice, but share this lattice with a dilutant without quadrupole moment: argon in the case of interacting 
N%, parahydrogen in the case of interaction with orthohydrogen, and KBr in the case of KCN , etc. Q. 

To construct the mean-field theory of a set of uniaxial quadrupoles interacting through randomly quenched and 
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frustrated isotropic exchange one may adopt the Hamiltonian 

H = — E E ^ItSlS^Sl — — E ^ 



(1) 



(i,3) 



where the spin vector S l is defined via the component /' = [SlSl — 8 pv /m) of the electric quadrupole moment 
tensor pH |. The summation runs over all the distinct pairs. Each S l has m components St (yU = 1, . . . , m) and, 
for convenience, is assumed to be a vector with fixed length |S l | = m. Of course, taking general m does not describe 
the experimental quadrupolar glasses. Rather it is a natural model to consider theoretically for classification. By 
analogy with SK, the spins are taken to interact via independent random interactions Jij which are assumed Gaussian 
distributed: 
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The mean J an d the variance J of the distribution depend on the total number of quadrupoles N to ensure a 
meaningful thermodynamic limit (N — ► oo) with an extensive energy: Jo = Jq/N and J = J/N 1 / 2 . The Hamiltonian 
JTJ can be seen as the Hamiltonian of an infinite-range model for N classical vector spins S* and zero external field 
is assumed. 

In Ref . [lfj it has been shown that in terms of the order parameters 
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tr [expL] -W>> 
the free energy per spin / - by means of the replica trick - is given by 
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where n is the number of replicas, [3 = l/fc^T, and L is 
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The elements of the order parameters Q°^ v \ p and M pv are not independent quantities and they can be parameterized 
in terms of five sets of independent parameters A a , B ab , C ab , D ab , and E ab . The non-zero extremal values of the 
above sets describe possible glass ordering [10). 

Upon decreasing the temperature, when one of the parameters becomes different from zero, the high-temperature 
disordered phase becomes unstable. By assuming continuous transition in the replica symmetric ansatz and provided 
that the average interaction is not too positive, 
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J 



< 



(6) 



there is a transition to an isotropic glass state occurring at the temperature Trs — (2mJ)/[fes(m + 2)]. This highest 
temperature phase transition is associated with the order parameter B ab acquiring a non-zero value, hinting at the 
presence of isotropic quadrupolar order |10|. As GS showed, below Trs the replica symmetric solution isotropic 
glass phase is unstable with respect to fluctuations in the space of broken replica symmetric isotropic glass order 
parameters, the instability being stronger than in conventional systems. 

It is worth noting that Eq. JBJ is the equivalent of the condition Jo/J< 1 in the SK model to ensure a transition from 
the paramagnet state to spin glass state. Furthermore, when the numerator of the right hand side of Eq. becomes 
negative, it is necessary to introduce a negative value for Jo in order to obtain a glass phase at low temperatures; this 
happens for m > 3.37 . . .. Increasing the negativity of J reduces the temperature of ferromagnetism onset, though it 
cannot be stopped. 
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III. ONE-STEP REPLICA SYMMETRY BREAKING THEORY 



The further analysis of GS leads to the conclusion that a replica symmetric ansatz cannot give a stable, and hence 
physical, solution of the quadrupolar glass model. Thus, one has to resort to a replica symmetry breaking ansatz. 
Close to the isotropic quadrupolar glass transition, i.e. confining attention to regions of parameter space (Jo, J) in 
which the highest transition temperature does correspond to a phase transition in the order parameter B ab , it is 
sufficient to consider only this parameter different from zero in the free energy (see Ref. |To|'). Thus, one requires Jo 
to satisfy the inequality JHJ for T < T#s, i.e. in the neighbourhood of the transition temperature. The free energy is 
then given by 



lim i [M!(m - + 2)£' (B^f + m (J3J) 2 £ V 

k ab ab 
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logTr exp 
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The symbol Yl'ab stands for a sum which excludes terms with any equal indices; i.e. a b. The paramagnetic 
contribution (3fpM — — (P J) 2 m(m — 3)/2 has been subtracted as it does not depend on the order parameter. One 
then looks for a replica symmetry breaking ansatz for B ab . Here we consider this to the first level in the standard 
Parisi procedure. 

Using the standard procedure of the replica symmetry breaking method, one groups the n replicas in blocks of 
x, where a; is a parameter (between 1 and n) to be located by the saddle points equations. Each block contains x 
replicas. Thus, one has 



B 



I'/ if I(a/x) = I(b/x) 
) otherwise 



(8) 



where I(y) is an integer valued function: its value is the smallest integer which is greater than or equal to y. Upon 
substituting in Eq. J7J, one has 
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where there are n/x blocks labelled by k — 1, . . . ,n/x. The index a belongs to block k if: I(a/x) — k. We shall 
now perturbatively compute the free energy @ by Taylor-expanding the free energy up to fourth order in q. This 
approximation implies that one is assuming the transitions to be continuous or at worst weakly discontinuous. After 
some lengthy but straightforward algebra, one obtains the following expression for the free energy: 
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where the coefficients in the free energy expansion are given by 

(to- 1)(to + 2) 3 
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4m 3 (1-tr . 
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and t is the reduced temperature, defined as t = 1 — (T 2 /Tq), where 

Tc = (JTW (12) 

For m less than a critical value to*, discussed below, the transition is continuous, with the same Tc as predicted by 
RS theory I n this case one may use the saddle point method to evaluate the extremal values of the parameters 
q and a. These saddle point equations are 

dq dx 

from which one has 

= 3ta 2 + 3q{a 3 + (x - 2)d 3 ) + q 2 [Pi + (a - 2)(/3 3 + (a: - 3)/3 5 )] 

= 6ta 2 + 4q{a 3 + (2a - 3)d 3 ) + <? 2 [Pi + (2a - 3)/3 3 + (3a; 2 - 12a + 11)/%] . ' 

According to the value of the quadrupole dimension to, the transition from the higher temperature RS phase can 
occur either continuously or discontinuously. The transition is continuous in q for to < m* ~ 3.37, i.e. when the 
coefficient 0:3 becomes larger than S3, q and x satisfy the following equations that express the extremum of the free 
energy functional ljTU|) 

t<0 : \ q = ° (14a) 
x undetermined 

'q = - 9 J a2 - < + Aa ,- t2a2 Vi -2 {{al + 4a 3 d 3 - 15d§) /? 5 

t>0 : <j + d 3 [-5/3id 3 + &(7d 3 -2a 3 )]} + 0[t 3 ] , (14b) 

_ 33 , to2[al(-ft+fa)+fe(^-2a 3 a3-a3)] , nr+2l 

~ &3 H 4(a 3 -d 3 )a^ h ^ J - 

Close to the transition temperature Tc, i.e. when t •> 0, the solution (|14b|l is valid only for 2 < to < m* ~ 3.37 within 
1RSB subspace. At to* the cubic term in the free energy functional (|10|l changes sign. This coincides with x = 1. 
Thus, above to*, the transition is not anymore continuous. The parameters q and x are plotted for m — 3 within this 
approximation in Fig. ^ By substituting Eq. <|14bl) in Eq. 1100 one finds the free energy of the glass phase close to 
the transition, to be 

24a 3 (a 3 — a 3 ) 96(m 2 — m — 8) 

Below the transition the free energy is larger than that of the paramagnetic phase. 

On the other hand, when to > to* there is still a glass solution to the saddle point equations but with a discontinuous 
onset of q from the higher temperature RS phase. The transition may be found with the additional requirement that 
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FIG. 1: (a) Plot of the one-step parameter q as a function of the temperature T (in units of fcs/J); (b) Plot of the one-step 
breaking parameter a; as a function of the temperature T (in units of fes/J). Both plots are drawn for m = 3 and within the 
quartic approximation for the free energy il l lit . 



the free energy in the paramagnetic phase is equal to the one in the glass phase with break point x equal to 1. 
Denoting this transition temperature by Td, to the quartic order in q for the free energy, one finds that 

t = 2(q 3 - a 3 ) 2 = 2(m + 6)(m 2 -m-8) 2 

D 3a 2 (A -f3 3 + 2/3 5 ) 9(m + 4)(168 + 34m-35m 2 -5m 3 +m 4 )' 1 ' 

where t D = 1 - (T D /T C ) . There is a discontinuous jump in the order parameter at Td from zero to 

2(q 3 -q 3 ) 2m(m + 6)(m 2 - m - 8)(1 - t D ) 

q ° /3i-/3 3 + 2/3 5 3(m + 2)(168 + 34m-35m 2 -5m 3 + m 4 )' 1 j 

where qjj = q(T = T^). In the neighbourhood of the transition temperature Tjj one finds that (x — 1) oc (t — tjj) and 
that the free energy of the quadrupolar glass phase is given by 



Pf = q D {t-t D f + 0[(t-t D ) 

Even though the transition is discontinuous in the order parameter q, there is no discontinuity in any thermodynamic 
quantities. Moreover, there is no latent heat at the transition. This behaviour is qualitatively common to a whole 
class of mean-field models of spin glasses e.g. the p > 2 spin model beneath a critical field Q and the Potts glass 
model above the critical Potts dimension p = 4 0, [a, ^| . 

Since the order parameter has a discontinuous jump at the transition temperature when m > m*, the perturbative 
approach should not be valid anymore. However, one may control the approximation by setting m = m* + e, where 
E<Cl. Thus, one obtains a quadrupolar glass phase with broken replica symmetry appearing below irj oc e 2 , with 
qo oc e and x(T — > Tp) — * 1. Explicitly, to leading order, one has 

2m*g(m*)(m* +4)(1 - 2m*) ^ r 9l , 

q D {m = m* + e) = -—^ e + 0[e% (18a) 

m* +2 

t D [m = m* + e) = 4.g(m*)(m* 2 - m* - 5/2) e 2 + 0[e 3 ], (18b) 

where g(m*) = (m* +6)/ [3(m* + 4)(m* 4 - 5m* 3 — 35m* 2 + 34m* + 168)] < 0. In the next section, we shall investi- 
gate the stability of the 1RSB solution found against small further RSB fluctuations. 

IV. STABILITY ANALYSIS 

In order to study the stability of the 1RSB ansatz one introduces li?5B-breaking fluctuations 

B ab ^ q S GaGb +f 1 ab , (19) 



\3 
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and expands the free energy to second order in the fluctuations n ab [l3|]. The group Kronecker delta So a G b is unity 
if a and b belong to the same group and zero otherwise [6j. One has to compute the second derivatives of the free 
energy JJJ) with respect to {B ab } at the 1RSB solution. The Taylor expansion of the free energy around the 1RSB 
solution is 



The quadratic form 



Bf f) 2 f 

f = f(qn S GaGb ) + v ab + dB J dBcd V ab V cd + ■■■ (20) 



9 f if b n cd (21) 



8B ab dB cd 



should be positive definite for a stable solution of the problem. It is easy to see that the only non-vanishing terms 
are the ones where a, 6, c and d all belong to the same group and the ones where a and c belong to group k, and b 
and d belong to group k'. The Hessian matrix S associated with the quadratic form 121|) factorizcs in n/x identical 
submatrices of dimension x(x — l)/2 x x(x — l)/2 which couple intragroup fluctuations and {n/x) [(n/x) — 1] /2 
identical submatrices of dimension x 2 x x 2 which couple intergroup fluctuations. We shall give details only on the 
intragroup matrices, because lengthy computations show that the intergroup matrices always have positive eigenvalues 
in the range of validity of the solution (see also Ref. @). 

If the 1RSB is stable thermodynamically, all of the eigenvalues of the stability matrix should be positive. The 
intragroup matrices S^ ab ^ cd ^ have three different types of matrix element: 

s (ab)(ab) = ^ 2 a 2 t - 4a 3 B ab - 2{3 1 (B ab ) 2 - -(3 3 ^ B bk B ka 

k=£aj=b 

s (ab)(ac) = _2a 3 B bc - -j3 3 [(B bc ) 2 + 2B ab B bc + 2B bc B ca ] 

- E Bk " BCk 

o 

g(ab)(cd) _ _Zp 5 (^B bc B da + B ac B db ) , 

3 

where a ^ b ^ c ^ d. Since a, 6, c and d belong to same group 

S (ab)(ab) = _ 2a2t _ m3q _ 2(3iq 2_^ x _ 2)q 2 
S (« 6 )(.c) = _ 2d 5 2_2 

6 3 

S (ab)(cd) = 

o 

The eigenvalues of the intragroup matrices, to order q 2 , are given by 



Xi = -2ta 2 - Aq [a 3 + (x - 2)a 3 ] - 2q 2 [0i + (x - 2)/3 3 + (x 2 - 5x + 6)/3 5 ] 

1 
6 



A 2 = -2ta 2 - 2q [2a 3 + (x - 4)d 3 ] - - q 2 [l2/3 x - (24 - 7x)(3 3 + 4(x 2 - $x + 18)/%] (22) 



A 3 = -2ta 2 + 4q(a 3 - a 3 ) - i q 2 [6ft + (a; - 7)/3 3 - i(x - 4)/3 5 ] . 

The behaviour of the eigenvalues in the ordered phase is obtained by substituting in to the above equations the values 
of the parameters q and x pertinent to the continuous or discontinuous transition. Close to the continuous transition 
temperature, one finds that the first two eigenvalues are positive in the range of validity of the solution 2 < m < m*: 



The last one, to order t 2 , 



(m-l)(m + 2) 3 , 9l 

(m + 2) 3 (m 3 - 3m 2 - 10m + 12) r „ 

4m z (m z — to — a) 



(to - l)(m + 2) 3 (2m 5 + 15m 4 - 8?ti 3 - 104m 2 - 32m + 96) ., 



(23) 



A - 3 = IB 2, , fiV 2 w >-t 2 + 0[t% (24) 

16to z (to + b)(m z — m — 8) z 



7 



is positive only for to > to?, — 2.46. Thus, one Imds a iKiit staDle mean-tield theory 
only in the range to?, < to < to*. This lower limit is the same as given in Ref. Q and obtained by means of a 
complementary calculation based on the full replica symmetry breaking (FRSB) ansatz near Tq with a perturbation 
treatment. For to > to* the behaviour of the eigenvalues in the ordered phase is obtained by setting x — 1 and by 
substituting q with the value qn obtained in Eq. (|18a|l . Upon inserting these values, one easily finds that all the 
fluctuations around the ordered 1RSB phase are finite and positive: 

g(m*) (m* - 1) (to* + 2) 3 (-IOto* 2 + IOto* + 3) 2 
Ai = A3 = -5 e > 

m 3 (25) 
g(m*) (m* - 4) (m* - 1) (to* + 2) J (2m* - 1) (2to* + 5) 

A2 = -5 e > 0. 

to z 

Thus, within our approximation, one finds a 1RSB stable mean-field theory with a discontinuous transition when 

to > TO*. 



V. DYNAMICAL TRANSITION 



Generally, disordered systems with a discontinuous transition have a temperature Tq where a dynamic instability 
appears. This temperature is called the glass temperature and is higher than the transition temperature To where 
the replica symmetry breaks thermodynamically, if the latter breaking is discontinuous. 

In the soft spin version of the Potts glass model [l4| it has been shown - by means of dynamical studies of the 
mean-field theory - that indeed there is another transition at temperature Tq > Td as in the p-spin model for p > 2 
[Tlj . Both static and dynamic transitions in the Potts (p > 4) case, have also been found in Refs. 0,0,01- I n the 
study of the thermodynamics of the quadrupolar glass, Tq can be computed by means of marginal stability ^||- By 
requiring the vanishing of the first and second derivative of the free energy 1)10(1 with respect to q 

y / q=qo v y / q=qc 

one finds, within our approximation, the dynamical transition temperature Tq and the corresponding discontinuous 
jump qo of the quadrupolar glass model 

= 3(q 3 - Q3) 2 = (to + 6)(to 2 -to-8) 2 

G Aa 2 (Pi -(3 3 + 2/3 5 ) 4(to + 4)(168 + 34to-35to 2 -5to 3 + to 4 ) 1 ' 

and 

3(q 3 - q 3 ) to(to + 6)(to 2 - to - 8)(1 - t G ) 

qG ~ 2(/3i - fa + 2(3 5 ) ~ 2(to + 2) (168 + 34to - 35m 2 - 5to 3 + to 4 ) ' [ ' 

where to = 1 — {Tq/Tc) and qc = q(T = Tq ). Again, by assuming the jump qo near the temperature Tq to be 
small, one can control the approximation by letting to = to* + e, e <C 1. Thus, one has 

! )(l-2m* 
2(to* +2) 



. „ . 3to* 5 (to*)(1 - 2to*)(to* +4) ^ 2l . . 

q G (m = m*+e) = — ^ — e + 0[e 2 ), (29a) 



t G (m = to* +e) = ^g(m*) (to* 2 - to* - 5/2) e 2 + 0[e\ (29b) 



Within the approximation used, one finds 



— = 7 + ^9(™*) (~2to* 2 + 2to* + 5) e 2 + O [e 3 ] . (30) 
qo 4 16 v 

It is worth noting that, to the leading order, exactly the same value for this ratio has been obtained for the Potts 
glass in Refs. 0,S^3' suggesting a sort of universality related to the same general structure of the free energy for 
both quadrupolar and Potts glass models (see also Ref. 0). The results for qo and q G as a function of l/log(m) are 
shown in Fig. [3 The ratio between the two transition temperatures Tq/Td is very close to one 

^ = 1 + 1 g{m*) (~2m* 2 + 2m* + 5) e 2 + O [e 3 ] , (31) 
Id 8 1 
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FIG. 2: The static qo = q(T = T D ) and dynamic qc = q(T = T G ) order parameter as a function of l/log(m). The dotted 
line is for the static value, the bold line is for the dynamical one. 



with Tq always bigger than To, as shown in Fig. For m — 4 (the smallest integer value of m compatible with a 
discontinuous transition) one finds that 

q D = 0.056... ; Td = 1.342... 

q G = 0.042... ; T G = 1.343... ^ ' 

Of course, a naive extension of our result to include large m results is not possible since we have assumed all the 
transitions to be at worst weakly discontinuous, implying the possibility to explore only the range m close to m*. 
However, since the mean-field theory of the Potts glass is qualitatively very similar to that of the quadrupolar glass, 
one may have an idea of the large m limit by considering the large p limit in the Potts glass model. There, the ratio 
Qg/qd stays close to 3/4 for a large range of p, though it increases for very large p, approaching unity. The ratio 
Tg/Td grows very slowly with p pj. This behaviour is different from the case of the p-spin model, where the ratio 
IgIid is not close to 3/4 even for small values of p and, in the limit p — > oo, it converges to unity. Moreover, the 
ratio Tq/Td grows faster with p than in the Potts problem 0. 

A useful and new representation of the phase diagram is obtained by plotting the phase boundary line in the plane 
(T, 1/m). One may identify three stable phases, RS paramagnetic and both 1RSB and FRSB glasses. In Fig. [31 
the phases are labelled by their symmetry breaking and the manner of the onset from the paramagnet. The 1RSB 
transition is continuous between m\ < m < m*, whereas it is discontinuous above m*. Note that, at m = m* , the 
transition from RS passes continuously from continuous 1RSB (C1RSB) to discontinuous IRSB (D1RSB) within 
the one-step RSB phase. The dotted line in the figure corresponds to the m > m* dynamical transition temperature 
given in Eq. (|27|l . The situation is analogous to that of the Potts glass model which shows a crossover from the 
continuous transition to the discontinuous transition as the number of Potts states increases [10. Thep>2 -spin 
Ising and spherical spin glasses also show transitions from C1RSB to D1RSB as an applied field is reduced but differ 
from the present problem in that the critical field makes also a maximum in the transition temperature, in contrast 
to the present monotonic variation with 1/m. For m < the transition is continuous to full replica symmetry 
breaking. A phase line (not shown, but continuous) separates the one-step and full replica symmetry breaking phases 
within the RSB region. 



VI. CONCLUDING REMARKS 



In this paper we have investigated the quadrupolar glass model in the framework of the replica method. Upon 
introducing a simple one-step replica symmetry breaking ansatz, one may find a stable mean-field theory with a 
continuous or discontinuous transition, according to the value of the quadrupole dimension m. The transition is 
continuous to one-step replica symmetry breaking in the range of the quadrupolar dimension 2.46 < m < 3.37. For 
the discontinuous transition (m > 3.37) there are two different transition temperatures. We have computed the ratio 
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1.4 



RS 




1/m 



FIG. 3: Phase diagram of the quadrupolar glass model in the (T, 1/m) plane. T is in units of fcs/J. For 2 < m < mf. the 
transition from the higher temperature RS phase occurs within the full replica symmetry breaking mechanism. For m < m* 
the thermodynamic and dynamical transitions coincide. For m > m* the dynamical transition, denoted by the dotted line, 
is higher that the thermodynamic one (solid line). The plot is shown within the quartic approximation for the free energy 
expansion in q. 



Qg/id, where qc and qd are the order parameters associated, respectively, with the dynamic and thermodynamic 
transition. The ratio between the two transition temperatures Tq/Tjj is also computed. Within the approximation 
used, the values of these ratios, qdlD = 3/4 and Tq/Td — 1 (to leading order), are the same as those found in Refs. 
[Ill III for the Potts glass model. 

The results we have obtained confirm the general wisdom that the properties of the quadrupolar glass, with con- 
tinuous (to < to*) and discontinuous transitions (to > to*), are similar to those of the p < 4 and p > 4 Potts glass 
well studied in the literature [a, S LDi ll3 ll3 • Although the investigation focused on the quadrupolar glass phase, in 
the wider (J , J,T) space there should exist different types of ferromagnet, collinear and canted, see e.g. Ref. [Toj . 

The full phase diagram should include also another curve which may be captured by complexity arguments. In 
analogy with the p-spin spin glass there should be another critical line T compL associated with the onset of macroscopic 
complexity [l7j . 
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